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Angular momentum can be transferred to a trapped Bose- 
Einstein condensate by distorting its shape with an external 
rotating field, provided the rotational frequency is larger than 
a critical frequency fixed by the energy and angular momen- 
tum of the excited states of the system. By using the Gross- 
Pitaevskii equation and sum rules, we explore the dependence 
of such a critical frequency on the multipolarity of the excita- 
tions and the asymmetry of the confining potential. We also 
discuss its possible relevance for vortex nucleation in rotating 
traps. 

PACS numbers: 03.65.-w, 05.30.jp, 32.80.-t, 67.40. Db 



It is well known that in order to create excitations in 
a superfluid by moving an impurity one has to overcome 
some critical velocity, which is fixed by the spectrum of 
the excited states of the system. In the case of a rectilin- 
ear motion of a heavy object, conservation of energy and 
momentum yields the critical velocity Vc = min[e(p)/p], 
where e(p) is the energy of an excitations carrying mo- 
mentum p ; this is known as Landau criterion lit]. The 
external perturbation can also be a moving boundary, 
Ukc the wall of a container, which produces excitations 
through its roughness. In liquid helium such a critical ve- 
locity has been the object of a longstanding investigation, 
involving the structure of the phonon-roton branch and 
the nucleation and dynamics of quantized vortex lines 
and rings. The existence of similar critical velocities in 
trapped Bose-Einstein condensates of alkali atoms is cur- 
rently under study. In Ref. |3] an external potential (a 
laser beam) is used to produce a hole in the conden- 
sate and this hole is moved back and forth at variable 
velocity, playing the role of a massive impurity. Theo- 
retical interpretations in terms of the Landau criterion 
have been already presented in Refs. ^^- Here we con- 
sider the analog process in which the perturbation acts 
at the surface of the condensate, producing a rotating 
shape deformation. This kind of perturbation has been 
recently implemented in Q to excite quadrupole and oc- 
topole modes, but the analysis is also relevant for rotat- 
ing traps, like the one used in Ref. [p| to produce quan- 
tized vortices. The similarity with liquid helium is again 
worth stressing; in that case, when the liquid fills a rotat- 
ing bucket, the roughness of the wall can transfer angu- 



lar momentum to the superfluid, favoring the occurrence 
of quantized vorticity above some critical angular veloc- 
ity H. However, since the crucial length scales (healing 
length, surface thickness, interatomic distance, etc.) in 
liquid helium and trapped condensates are different, the 
role of surface excitations on the Landau criterion is also 
very different. In particular, surface modes can be ig- 
nored in the analysis of the experiments with helium in 
rotating vessels, while they might be crucial for gases of 
alkali atoms in rotating traps, as discussed below. 

Let us consider a condensate which is initially in its 
ground state in an axially symmetric trap, with trapping 
frequencies uj± and uj^ = XoJ± , where A is the anisotropy 
parameter of the harmonic confining potential. Then, let 
us imagine a small distortion of the condensate caused 
by an external perturbation which rotates around the z- 
axis at angular velocity f2. This perturbation may be 
a modulation of the confining potential, a rotating laser 
beam, or something else; we only assume here that the 
consequent distortion of the condensate is indeed a shape 
deformation, coupled to the surface modes of the system 
(ripplon-like waves) . Let us consider the process in which 
the external probe produces a single excited state, cor- 
responding to a surface wave with energy hLo{£) and an- 
gular momentum £ along z (azimuthal quantum number 
m — £). One can assign a moment of inertia / to the ro- 
tating perturbation and write the conservation laws for 
energy and angular momentum, letting / ^ 00 at the 
end. The simple result is that the conservation laws are 
both satisfied only if fi is larger than the critical an- 
gular velocity flc — nimluj {£)/£]. This is the analog of 
the Landau criterion. If the rotation is slower than flc, 
the condensate simply adjusts its shape to the rotating 
perturbation by moving as an irrotational and nondissi- 
pative fluid. Conversely, if the rotation is faster than flc, 
the external perturbation can excite surface modes in a 
continuous way, transferring energy and angular momen- 
tum to the condensate. 

The excitation spectrum of a condensate can be 
accurately evaluated by solving the linearized Gross- 
Pitaevskii (GP) equation, also known as Bogoliubov's 
equations |q| . The theory provides the energy of the ex- 
cited states at zero temperature in terms of the relevant 
quantum numbers (number of radial nodes and angu- 
lar momentum). For our purposes what matters is the 



energy hu!{£) of the surface modes, i.e., the ones with 
no radial nodes. In a previous work, we already pre- 
sented some results for the critical angular velocity flc 
in a spherical condensate B. Here we discuss more sys- 
tematically the consequences of this critical behavior and 
provide results for nonspherical traps. 

We avoid the complication of solving the Bogoli- 
ubov's equations for axially symmetric condensates by 
using sum rules. As already shown in pO| , the sum 
rule approach provides rigorous upper bounds to the 
dispersion uj(£) using the ground state density, solu- 
tion of the stationary GP equation, as main ingredient. 
The key quantities are the fc-energy weighted moments, 
TOfc = /„ SF{E)E^dE of the dynamic structure fac- 
tor, Sf[E) = J2j \{J\P\Q)?5{E-Eja), associated with a 
given excitation operator F , where Ejq — (Ej — Eq) is the 
excitation energy of the eigenstate \j) of the Hamiltonian. 
The TOi and m^ moments can be written as expecta- 
tion values of commutators, mi = (l/2)(0|[Ft, [H,F]]\0) 
and mg = (l/2)(0|[[Ft, i/], [iJ, [iJ,F]]]|0), and the ratio 
(ma/mi)^/^ gives the desired upper bound to the energy 
of the states excited by the operator F. In order to select 
the surface modes, one chooses the excitation operator in 
the form F ~ J2i=i ^lYu{9i,<j)i), with £ 7^ 0, where N is 
the number of atoms in the condensate. The explicit cal- 
culation of the ratio [m^/mi) with this operator gives 
the dispersion law 



uj^{e)^ujl£[l + {£^l)/3t^] 
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and rj_ = {x^ -\- y^)^/^. Here n{r±_^ z) is the ground state 
density of the condensate which can be calculated by 
solving the stationary GP equation. The sum rule es- 
timate is close to the exact Bogoliubov's spectrum when- 
ever the state excited by F is highly collective, with a 
strength which almost exhausts the dynamic structure 
factor. This is true for shape deformations of low mul- 
tipolarity, which are the important ones for the present 
analysis. The dynamic structure factor associated with 
large £ excitations is expected to be more and more frag- 
mented and, consequently, the sum rule upper bound be- 
comes less accurate. 

The surface state with lowest energy is the quadrupole 
mode. The sum rule estimate is easily obtained by in- 
serting £ — 2 into Eqs. (0)-(||). One finds 
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where (£'kin)_L and (£'ho)_L sue the radial contributions to 
the kinetic and external potential energies, respectively. 
For noninteracting particles these two energies are equal 



and one gets w = 2wj_. The other interesting limit is 
N —f 00 (Thomas- Fermi limit), when the kinetic energy 
in the ground state is much smaller than both the atom- 
atom interaction energy and the trapping potential. In 
this limit Eq. (^) gives to — \f2 w_l. The latter result can 
also be obtained by solving the hydrodynamic equations 
of a supcrfluid |TQ] , which provide the general dispersion 
law w(£) = vt ijj}_ for the surface modes with m = ±^. 
The hydrodynamic approach, however, is accurate only 
for the lowest multipoles when one can completely ne- 
glect the contribution of the quantum pressure term in 
the GP equation; for the values of N which are rele- 
vant in the current experiments, the hydrodynamic dis- 
persion becomes soon inaccurate as £ increases. This has 
an important consequence for the evaluation of the criti- 
cal frequency f2c — min[ti> (£)/£], since the hydrodynamic 
dispersion gives uj{€)l£ ~ uj^/\/£ and hence Vic — 0; con- 
versely, Vic always exhibits a minimum, even for very 
large iV, when calculated with Bogoliubov's equations 
or using the sum rule results (|l|)-(g). 

In Fig. H we show typical results for uj{£)/£, in units 
of uj^, obtained by inserting in Eqs. (|l|)-(§) the ground 
state density of three different condensates of ^^Rb atoms 
(we use a = lOOao for the scattering length, where ag is 
the Bohr radius). The stationary GP equation has been 
solved numerically as in Ref . |ll| . The condensates have 
the same number of atoms, N — lO'^, and the same aver- 
age trapping frequency, Who = {^x^yi^JzY'^ = \^'^0Ji_ = 
27r(20Hz), but different asymmetry parameter: A = 0.1 
(prolate), A = 1 (spherical) and A = 10 (oblate). The iV- 
dependence of bj{£)/ £ is instead given in Fig. 0, where we 
plot the curves for condensates with different N but the 
same prolate geometry: A = 0.06 and Who = 27r(82.5Hz), 
close to the values of Ref. 0. In order to test the ac- 
curacy of the sum rule approach, we have also solved 
numerically the Bogoliubov's equations for a spherical 
condensate. The results are shown as open circles in 
Fig. [l| and can be compared with the sum rule results 
for A = 1 (solid circles). As one can see, the accuracy of 
the sum rule approach is very good for the lowest values 
of £, remaining of the order of 10% up to ^ = 10. This 
is enough for our purposes. Finally, the hydrodynamic 
prediction is also shown as a dashed line in both figures. 

In Table I we summarize the situation for three dif- 
ferent condensates. The parameters of the first (prolate) 
and the second (spherical) are close to the ones of current 
experiments at ENS-Paris and JILA-Boulder [n2[, re- 
spectively. The third one (oblate) is presented in order 
to stress the dependence of the critical frequencies on the 
geometry of the trap. Actually, the critical frequencies 
ijj{£)/£ for the lowest multipoles {£ = 2,3,4) turn out 
to be rather independent of the shape of the condensate 
and are well described by the hydrodynamic dispersion, 
which is independent of A. For the quadrupole mode, 
for example, it is very well reproduced by the hydrody- 
namic value Lo±/^/2 = O.TOTwj,. The minimum oiijj{£)/£ 



instead depends on A in a more significant way, since it 
involves higher ts, especially for the oblate geometry. 

The frequencies uj{£)/£ in Table I are thresholds for the 
transfer of angular momentum to a condensate through 
shape deformations. The values of £ actually involved in 
a given experiment depend on the type of process used. 
For instance, a laser beam which makes a rotating hollow 
on the surface of the condensate might excite a bunch of 
states with several £'s, just as happens for a hole mov- 
ing through the fluid and exciting phonon-like excita- 
tions with several momentum components. Conversely, 
the stirred confining potential of Ref. B, with a rotat- 
ing quadrupole deformation, is expected to excite mostly 
quadrupole excitations. What really happens just above, 
or close to the threshold can not be described only in 
terms of elementary excitations. In fact, once a thresh- 
old is reached, the rotating perturbation can pump many 
excitations into the condensate, bringing it far from the 
linear regime. A wider configuration space is made avail- 
able and the system can find new rotating equilibrium 
states, like the stationary deformed configurations dis- 
cussed in p3| , or jump into states with quantized vortic- 
ity. 

Concerning vortices, an important frequency scale is 
fixed by the difference between the energy per particle 
of a condensate with and without vortex calculated in 
the rotating frame, fl^ = {l/h)[{E/N)y ~ {E/N)g\. This 
corresponds to the ratio between the 'excitation' energy 
of the vortex in the rest frame, {Ey —Eg), and its angular 
momentum Nfi. When the rotation frequency fl is equal 
to ^y the energies of the configurations with and without 
vortex are equal in the rotating frame and, above il^,, 
the vortex state becomes the global energy minimum. 
In this sense, fly is analog to lu{£)/£, since the latter 
represents the frequency above which a surface state of 
angular momentum £ has negative energy in the rotating 
frame |14f| . A crucial difference, however, is that in order 
to transfer Nh angular momentum and nucleate a vortex 
above Qy one has to overcome an energy barrier, while no 
barrier is found above uj{£)/£ for the creation of quanta 
of angular momentum through surface excitations. 

The quantity fly can be obtained from the solution of 
the stationary Gross-Pitaevskii equation for the ground 
state and the state with a vortex line along z. We have 
already presented some numerical results for this quan- 
tity in ||ll[. By repeating the same kind of calculation 
for the three condensates in Table I we find the values 
given in the last column |13] . It has been noted in Ref. ra] 
that the experimental result for the frequency at which 
vortices occur is significantly larger than Qy. The crit- 
ical frequency in the experiment turns out to be about 
0.67a;j_, to be compared to the value 0.35a;_L given in 
the first row of Table I. This is likely to be the con- 
sequence of the above mentioned energy barrier, which 
prevents the nucleation of vortices even above fly . If one 
tries to impart a rotation to the condensate by rotating 



a slightly asymmetric external field, the crucial question 
is how angular momentum is actually transferred to the 
system. Our analysis suggests that, if the transfer mech- 
anism is the creation of surface excitations of multipo- 
larity £, the relevant threshold is uj{£)/£, which can be 
significantly larger than Qy, as shown in Table I. This 
argument might be used to interpret the fact that the 
experimental critical frequency of Ref. turns out to 
be closer to the quadrupole value w(2)/2 ^ uj^/^/2 than 
to Uy. Were this true, one should find almost the same 
critical frequency by repeating the experiment in spher- 
ical or oblate condensates. This clearly suggest that the 
A-dependence could discriminate between alternative de- 
scriptions for vortex nucleation and stability (see, for in- 
stance, Refs. |l6|-[l8| and references therein). Another 
important signature of the suggested threshold mecha- 
nism would be the observation of strongly deformed con- 
figurations in the initial stage of the vortex nucleation. 
Preliminary evidence for such deformations has already 
emerged ||l^. Finally, the threshold behavior here dis- 
cussed would imply a significant hysteresis when the ro- 
tation frequency is first speeded up to produce vortices 
and then slowed down again, since one expects differ- 
ent mechanisms for the two processes: the nucleation of 
the vortex state at relatively high rotation frequency and 
the destabilization of the vortex at lower frequency (for 
an extensive discussion about the destabilization of vor- 
tices see [M). The hysteresis mechanism of the vortical 
configurations has some similarity with the behavior of 
superfluid helium. In that case, however, vortex nucle- 
ation is associated with the pinning of vortex lines at the 
walls of the container R] . Trapped condensates have no 
rough surfaces or, in other words, they have an adjustable 
"roughness" ; in this sense, they offer a novel and com- 
plementary tool for understanding nucleation processes 
of quantized vorticity in superfluids. 
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FIG. 1. The quantity ijj{l)/l, in units of lli±, obtained 
with Eqs. (llh-M as a function of £ for three condensates 
of *^Rb atoms with different values of A = loz/ujj_: pro- 
late (A = 0.1), spherical (A = 1), oblate (A = 10). The 



three condensates have the same value of A'^ 



10'' and 



i.^ho ~ \^ ujj_ — 27r(20Hz). Open circles correspond to the 
exact solution of the Bogoliubov's equation in the spherical 
case (A = 1). The dashed line is the hydrodynamic prediction 

uj{e)/e = i/VI. 
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prolate 
A ^ 0.0058 


0.72 


0.61 
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0.53 


0.35 


i'± = 175Hz 
N = 2.5x 10^ 












spherical 
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TV = 3 X 10^ 


0.71 
(0.71) 
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(0.59) 


0.53 
(0.52) 
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(0.41) 


0.24 


oblate 
A = 10 


0.71 


0.58 


0.50 


0.33 


0.12 


v^ = lOHz 
A^ = 3 X 10^ 













Table I. Critical frequencies for shape deformations in 
different condensates. First column: type of condensate. 
Next four columns: critical frequency uj{£)/£, in units of 
the radial trapping frequency U!± = 2ttv± , obtained with 
the sum rule expressions (|l|)-(§) for £ = 2, 3, 4 and for its 
minimum value. The numbers in brackets are the exact 
results of Bogoliubov's equations in the spherical case. 
The last column is the quantity fly — {l/h)[{E/N)y — 
{E/N)g\ (see text). 
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FIG. 2. Same as in Fig. 1 but for condensates with 
different values of A*'. All condensates have the same 
ojho = 27r(82.5Hz) and A = 0.06. The parameters used for 
the curve N = 2 x 10"" are close to the actual situation in the 
experiments of Ref. m. 



